BUNDfSREPUBLIK DEUISCH 



9 




bzy o 
V3 r ~ 



Oft 




03 DEC 2004 



Prioritatsbescheinigung iiber die Einreichung 
einer Patentanmeldung 




Aktenzeichen: 
Anmeldetag: 
Anmelder/lnhaber: 
Bezeichnung: 



102 24 735.8 



04. Juni 2002 




Dr. Christof H o I b e r g , Iffeldorf/DE 



Verfahren, Vorrichtung und Computerprogrammpro- 
dukt zur Erzeugung eines dreidimensionalen Modells 



IPC: 



G 06 T, G 06 F 




Die angehefteten Stiicke sind eine richtige und genaue Wiedergabe der ur- 
sprunglichen Unterlagen dieser Patentanmeldung. 



Munchen, den 11. Juni 2003 
Deutsches Patent- und Markenamt 
Der President 

Im Auftrag 



Horn 




A 9161 

03/00 

EDV-L 



CERTIFIED COPY OF 
PRIORITY DOCUMENT 



BEST AVAILABLE COPY 




Tft)TKE - BUHLING - KESH#fe PARTNER (GbR) 

TBK-Patent POB20 19 18 80019 MUnchen ~\ Patentanwalte 

Dipl.-lng. HarroTtedtke 
Dipl.-lng. Relnhard Wnne 
DiplMng. Hans-Bemd Pellmann 
DipMng. Klaus Grams 
Dipl.-lng. Aurel Vollnhals 
Dipl.-lng. Thomas J A. Leson 
Dipl.-lng. Dr. Georgi Chivarov 
DIpl.-lng. Matthias Grill 
Dipl.-lng. Alexander KQhn 
Dipl.-lng. Ralner Bdckelen 
Dipl.-lng. Stefan KJingele 
Dlpl.-Chem. Stefan BOhiing 
DlpWng. Ronald Roth 



4. Juni 2002 
DE 33107 



Dr. Christof Holberg 

82393 Iffeldorf, Deutschland 




VERFAHREN, VORRI CHTUNG UND COMPUTERPROGRAMMPRODUKT ZUR 
ERZEUGUNG EINES DRE I D IMENS I ONALEN MODELLS 



Dresdner Bank MOnchen Kto. 3939 844 BLZ700 800 00 
Deutsche Bank MOnchen Wo. 286 1060 BLZ 700 700 10 
Poslbank MOnchen Wo. 67043 804 BLZ 700 100 80 

DaWchWfcngyoBank DOsseldorf Kto. 8104233007 BLZ 300 207 00 
Sanwa Bank DOsseldorf Kto. 500 047 BLZ 301 307 00 



Telefon: +4989 544690 

Telefax (G3): +49 89 532611 
Telefax (G3+G4): +4989 5329095 
E-Mail: postofflce@tbk-patentde 
Internet http://toww.fok-patent.de 
Bavariaring 4-6, 80336 MUnchen 



BE SCHRE IBUNG 



Die Erfindung betrifft ein Verfahren, eine Vorrichtung 
und ein Computerprogrammprodukt zur Erzeugung eines drei- 
dimensionalen Modells fUr ein real bestehendes Objekt, 
insbesondere zur Erzeugung eines Flachen- oder Volumen- 
5 korpermodells oder eines FE-Modells (FE: Finite Elemente) 
aus digitalisierten Daten des Objekts. 

Auf vielen technischen Gebieten besteht der dringende 
Bedarf, das Abbild eines real bestehenden Objekts in ein 
10 moglichst wirklichkeitsnahes Computermodell zu uber- 

ftlhren. Wahrend sich einfache Objekte meist noch recht . 
gut mit Hilfe von CAD-Programmen (CAD: computerunter- 
sttttztes Konstruieren) , FE-Programmen oder dergleichen 
nachkonstruieren lassen, wird dies bei komplizierten 
15 Strukturen immer schwieriger. 

Besonders auffal.lig ist die Diskrepanz zwischen Bedarf 
und tatsachlich Machbarem in den Biowissenschaften, wenn 
beispielsweise die hochkomplizierte Gestalt von 

20 anatomischen Strukturen modelliert werden soil. Die 

geometrischen bzw. morphologischen Ungenauigkeiten des 
Modells kommen im besonderen Mafie dann zum Tragen, wenn 
auf der Grundlage des Modells eine FE-Analyse durch- 
geftihrt werden soil, um das physikalische Verhalten des 

25 Objekts zu berechnen. 

Bei der Erzeugung eines aus Knoten und Elementen 
bestehenden FE-Modells kann prinzipiell zwischen zwei 
Methoden unterschieden werden, und zwar zwischen der 
30 direkten und der indirekten Erzeugung des FE-Modells. Bei 
der direkten Methode werden die Knotenpunkte dem FE- 
Programm fest vorgegeben, wahrend sich das FE-Programm 
bei der indirekten Methode die Knotenpunkte aus vorge- 
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gebenen geometrischen Element en (z.B. aus Flachen- oder 
Volumenkorpern, Linien oder Punkten) selbst auswSLhlt. 

Untersucht man die bislang in den Biowissenschaf ten 
5 erstellten FE-Modelle, so fallt auf, dass die meisten der 
vorgestellten Verfahren direkt e Verfahren sind, bei denen 
dem FE-Programm die Knotenpunkte fest vorgegeben werden. 

Fig* 13 zeigt ein Beispiel ftir ein direkt durch 
10 Konstruktion erzeugtes FE -Mode 11 der Gesichtsweichteile, 
wie es M. Motoyoshi et al. in "Finite element model of 
facial soft tissue. Effects of thickness and stiffness on 
changes following simulation of orthognatic surgery", 
J Nihon Univ Sch Dent 35, Seiten 118-123 (1993) 
15 vorgestellt haben. Bei diesem Verfahren wird auf eine 
morphologisch exakte Ubertragung der Objektstruktur in 
den virtuellen Raum verzichtet und es wird versucht, die 
Komplexitat der Struktur so gut wie mSglich durch 
manuelle Nachkonstruktion nachzuahmen. 

20 

FE-Modelle konnen auch direkt durch schichtweise 
Vernetzung erzeugt werden. Dabei werden die geometrischen 
Daten des Objekts durch ein schichtdiagnostisches 
Verfahren oder durch Anfertigen histologischer Schnitte 
25 gewonnen, 

Bei der punktbasierten, schichtweisen Vernetzung werden 
die gewonnen Schichtbilder digitalisiert und die Grenzen 
der interessierenden Struktur identif iziert . Anschliefiend 

30 werden in jeder Schicht auf den Grenzlinien Knotenpunkte 
definiert/ die zuerst zweidimensional in der jeweiligen 
Schichtlage und anschliefiend dreidimensional zwischen den 
einzelnen Schichtlagen vernetzt werden. Fig. 14 zeigt ein 
solches durch punktbasierte, schichtweise Vernetzung 

35 erzeugtes FE -Model 1 eines Zahns, wie es C. Lin et al. in 
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"Automatic finite element mesh generation for maxillary 
second premolar", Comput Methods Programs Biomed 59, 
Seiten 187-195 (1994) beschrieben haben. 

Die schichtweise Vernetzung kann auch voxelbasiert 
erfolgen. Dazu wird tiber jedes gewonnene Schichtbild ein 
definiertes quadratisches Gitternetz gelegt und jedem 
Quadrat, das in der Schicht einem Voxel entspricht, ein 
kubisches Element zugeordnet. Jedes Element, das in der 
darunter liegenden Schicht nicht zu einem gewissen 
Prozentsatz von der interessierenden Struktur bedeckt 
ist, fallt weg. Obrig bleiben Schichten aus gleich- 
formigen Elementen, die ubereinander gestapelt ein 
dreidimensionales FE -Mode 11 aus kubischen Elementen 
bilden. Fig. 15 zeigt ein solches durch voxelbasierte, 
schichtweise Vernetzung entstandenes FE -Model 1 des 
menschlichen Schadels, wie es D. Camacho et al. in "An 
improved method for finite element mesh generation of 
geometrically complex structures with application to the 
skullbase", J Biomech 30, 1067-1070 (1997) beschrieben 
haben. 

Einen etwas anderen Ansatz haben P. van Zyl et al. in 
"Three-Dimensional finite element of a human mandible 
incorporating six osseointegrated implants for Stress 
Analysis of mandibular cantilever prostheses", Int J Oral 
Maxillofac Implants 10, Seiten 51-57 (1992) gewahlt. Das 
in Fig. 16 gezeigte FE -Model 1 eines Unterkiefers wurde 
direkt durch dreidimensionale Vernetzung einer aus 
Ref lexmikroskopaufnahmen gewonnenen Punktwolke erzeugt. 
Die Punktwolke wird direkt in das FE-Programm eingegeben, 
wobei die Punkte der Punktwolke als Knotenpunkte fur das 
FE -Mode 11 herangezogen werden. Allerdings ist die 
automat ische Vernetzung der Knotenpunkte nicht besonders 
zuverlassig. Die* gangigen FE-Programme haben haufig 



Schwierigkeiten mit der Vernetzung unterschiedlich 
dichter Punktwolken, sodass in der Regel spezielle 
Programme zur dreidimensionalen Vernetzung der Punktwolke 
verwendet werden mtlssen. Da diese Programme nicht 
kompatibel zu den FE-Standardprogrammen sind, entstehen 
auf der Sof twareseite Insellosungen, die nur von einem 
Fachmann bedient werden kennen. 

Indirekte Verfahren zur Erzeugung eines FE-Modells haben 
sich in erster Linie in den Ingenieurwissenschaften 
durchgesetzt. Bei den indirekten Verfahren werden dent FE- 
Programm beliebige geometrische Elemente vorgegeben, aus 
denen sich das FE-Programm selbsttatig die Lage der 
Knotenpunkte errechnet, wobei die geometrischen Elemente 
lediglich die Rander und Grenzf lachen des spateren FE- 
Modells festlegen. Die geometrischen Elemente konnen 
entweder im FE-Programm selbst erstellt werden oder uber 
eine sogenannte CAD/FEM-Kopplung als Flachen- oder 
Volumenkorper aus einem CAD-Programm importiert werden. 

Urn die Koordinaten der Knotenpunkte zu erstellen, gibt es 
bei der indirekten Erzeugung des FE-Modells zwei 
MSglichkeiten: Beim sogenannten "Mapped Meshing" werden 
die Knotenpunkte so festgelegt, dass viereckige oder 
hexaederfSrmige Elemente gebildet werden. Im Gegensatz 
dazu werden beim "Free Meshing" dreieckige oder 
tetraederfSrmige Elemente mit Zwischenknoten (sogenannte 
parabolische Elemente) gebildet, die sich besonders gut 
an komplexe Geometrien anpassen. 

Der entscheidende Punkt beim indirekten Verfahren ist, 
dass die dem FE-Programm vorgegeben geometrischen 
Elemente in der Regel manuell im FE-Programm oder in 
einem CAD-Programm konstruiert werden mtissen. In den 
Ingenieurwissenschaften stellt dies keineh Nachteil dar, 
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da die meisten Produkte sowieso mit Hilfe von CAD- 
Programmen entworfen werden. In den Biowissenschaf ten 
allerdings hat sich das indirekte Verfahren bis auf 
Ausnahmen, etwa bei der Beurteilung konstruierbarer 
5 Fremdkorper wie Hiif tprothesen usw., nicht durchsetzen 
k6nnen. Die Ungenauigkeiten, die entstehen, wenn ein 
reales Objekt durch direkte Konstruktion nachgeahmt wird, 
sind einfach zu grofi. 

10 Den bisher einzigen Losungsweg, das Abbild eines real 
bestehenden Objekts in ein der CAD/ FEM-Kopp lung 
zugangliches Format umzuwandeln, bietet das Reverse 
Engineering. Beim Reverse Engineering wird das Objekt 
digitalisiert und werden mittels Flachenrekonstruktion 

15 CAD-Flachen erzeugt, die sich dem Grundsatz nach in ein 
FE-Programm importieren lassen. Die CAD-Flachen setzen 
sich zumeist aus frei fonabaren Bezier- oder NURBS- 
Patches (NURBS: Nicht-uni forme rationale B-Splines) 
zusammen, die iiber ein Netz aus Kontrollpunkten 

20 sttickweise an die Oberf lachenf orm des Objekts angepasst 
werden. Bei den NURBS-Patches handelt es sich in der 
Regel urn mindestens bikubische parametrische FlSchen- 
elemente, die die Objektflache jeweils durch zwei poly- 
nomiale Kurven dritten Grades approximieren* Das Losen 

25 der linearen Gl e i Chungs syst erne wird mit zunehmendem Grad 
der Polynome immer aufwendiger, erlaubt aber eine 
genauere Anpassung der Flachenelemente an die Ober- 
fiachengestalt des Objekts. Aufgrund der komplizierten 
Algorithmen, die bei der Flachenrekonstruktion zum 

30 Einsatz kommen, sind die zum Reverse Engineering 
verwendeten Sof twareprogramme bislang sehr teuer. 
Aufierdem sind die Berechnungen mit hohem Zeitaufwand 
verbunden. 
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Daneben hat sich herausgestellt, dass beim Reverse 
Engineering die Anpassung der B6zier- oder NURBS-Patches 
an die Oberf lachenform einer geometrisch komplexen 
Struktur oft sehr ungenau ist und dass dabei haufig 
5 irregulare CAD-Flachen erzeugt werden, die sich weder mit 
CAD-Programmen noch mit FE-Programmen weiterverarbeiten 
lassen. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren, 
10 eine Vorrichtung und ein Sof twareprogrammprodukt zur 

Verftigung zu stellen, mit denen sich mit relativ geringem 
Rechenaufwand ein genaues dreidimensionales Modell eines 
real bestehenden Objekts erzeugen lasst. 

15 Diese Aufgabe wird durch ein Verfahren gemafi Anspruch 1, 
eine Vorrichtung gemafi Anspruch 6 und ein Computer- 
programmprodukt gemafi Anspruch 9 gelost. 

Die Erfindung zeichnet sich dadurch aus, dass das 
20 interessierende Objekt zunachst digitalisiert wird, urn 
ein Netzmodell des Objekts zu erzeugen, das Netzmodell 
anschliefiend in bilineare Flachenelemente zerlegt wird 
und die bilinearen Flachenelemente schliefilich zu einem 
Fl&chen- oder Volumenk6rpermodell vereinigt werden. 

25 

Unter dem angesprochenen Netzmodell sind Vielf lichen-, 
Oberflachen- oder Polygonnetze zu verstehen, die 
typischerweise aus einer Menge endlich vieler Polygone 
bestehen, bei denen jeweils zwei Eck- oder Knotenpiinkte 
30 eine Kante definieren und mehrere solcher Kanten einen 
geometrischen Korper beschreiben. Die geometrische 
Beschreibung des Korpers erfolgt in dem Netzmodell rein 
numerisch, d.h. im Gegensatz zu einem analytischen Ansatz 
wird die geometrische Form nicht durch mathematische 
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Gleichungen definiert, sondern allein durch die Lage und 
Dichte der Eck- bzw. Knotenpunkte. 

Da solche Netze rein numerisch beschrieben werden, lassen 
sie sich durch Digitalisieren des Objekts erzeugen. Das 
Digitalisieren kann dabei auf unterschiedliche Weise 
erfolgen. So kann das Objekt beispielsweise optisch oder 
beruhrend abgetastet werden # so dass sich eine die 
Objektoberflache beschreibende Punktwolke ergibt. Aus 
dieser Punktwolke werden dann die Knotenpunkte fur die 
Netzdarstellung gewonnen. Von dem zu digitalisierenden 
Objekt kennen aber auch Oberflachen- oder Schnitt- 
aufnahmen angefertigt werden. Mit Hilfe dieser Aufnahmen 
werden dann die Grenzen (Boundaries) des Objekts 
identifiziert und werden einzelne Punkte dieser Grenzen 
wiederum als Knotenpunkte fur die Netzdarstellung 
verwendet . 



20 



Bislang ging man davon aus, dass die numerischen Daten 
eines solchen Netzmodells nicht ohne Weiteres dazu. heran- 
gezogen werden k6nnen, einen Flachen- oder Volumenkorper 
zu erzeugen, der durch analytische Funktionen beschrieben 
wird. So wird in der Literatur zum Beispiel davor 
gewarnt, dass in CAD-Programmen zwar Volumenkorper in 
25 Oberflachennetze konvertiert werden konnen, diese aber 
umgekehrt nicht in Volumenkdrper konvertiert werden 
kSnnen (vgl. Benutzerhandbuch AutoCAD 2000™, s. 719). Um 
numerische Daten in analytische Daten umzuwandeln, musste 
bxslang der aufwendige und insbesondere bei komplexen 
Geometrien nicht immer erfolgreiche Weg des Reverse 
Engineering eingeschlagen werden. 

Die Erfindung schlagt eine BrUcke zwischen numerischer 
und analytischer Beschreibung der Objektdaten, indem sie 
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die numerischen Daten des Netzmodells in die analytischen 
Daten bilinearer Flachenelemente zerlegt. 

Unter bilinearen Flachenelement en sind Fiachenstticke zu 
5 verstehen, die jeweils durch zwei polynomiale Kurven 
ersten Grades bzw. durch zwei Strecken definiert sind. 
Die Endpunkte der Strecken ergeben sich dabei aus den 
Knotenpunkten des Netzes. Die beiden Strecken jedes 
Fiachenelements bilden zwei Kanten eines Polygonzuges, 
10 dessen tibrige Kanten sich durch Verbinden der Strecken- 

endpunkte ergeben. Jedes Flachenelement hat seine eigenen 
Kanten, die es nicht mit den angrenzenden Flachen- 
elementen teilt. 

15 Die bilinearen Flachenelemente sind vorzugsweise 

dreieckig, da sich Dreiecksf lachen besonders gut an 
komplexe Geometrien anpassen lassen. Bei einer solchen 
Dreiecksf lache fallt daher jeweils ein Endpunkt der 
beiden Strecken des Fiachenelements in einem Punkt 

20 zusammen und wird die dritte Kante des Dreiecks durch die 
verbleibenden Endpunkte der Strecken gebildet. Das 
bilineare Flachenelement kann jedoch prinzipiell auch 
eine Form mit vier Kanten (z.B. ein Quadrat) einnehmen, 
bei der die Endpunkte der Strecken nicht in einem Punkt 

25 zusammenf alien. 

Sofern die bilinearen Flachenelemente durch ein CAD- 
Programm bearbeitet werden, sollten sie vorzugsweise in 
Form von NURBS-Patches vorliegen, da NURBS-Patches gegen- 
30 tiber Rotations-, Skalierungs-, Translations- und 

Projektionsoperationen invariant sind. Dabei ist zu 
beachten, dass die Moglichkeit, NURBS-Patches hoheren 
Grades frei formen zu kSnnen, nicht genutzt wird, da die 
verwendeten Flachenelemente nur bilinear sind. 



35 
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Um das Netzmodell in bilineare NURBS-Patches umzuwandeln, 
wird das Netz vorzugsweise in das IGES-Format (IGES: 
Initial Graphics Exchange Specification) konvertiert. Das 
IGES-Format ist ein ANSI-Standard, der ein neutrales 
5 Format fur den Datenaustausch zwischen unterschiedlichen 
CAD-, CAM- (CAM: rechnergestutzte Fertigung) und 
Computervisualisierungssystemen definiert. Die bilinearen 
Flachenelemente entsprechen dabei IGES-Elementen der 
Nummer 128, die ftlr rationale B-Spline-Flachen vorgesehen 
10 sind. 



15 
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Die einzelnen bilinearen Flachenelemente, in die das Netz 
zerlegt wurde, werden wieder vereinigt, um einen 
geschlossenen Flachenverbund oder einen geschlossen 
Volumenkorper zu erzeugen. Dies geschieht dadurch, dass 
die gegenuberliegenden Kanten von zwei benachbarten 
Flachenelementen zusammengeheftet werden. Es werden also 
die zuvor getrennten Kanten der Flachenelemente zu einer 
gemeinsamen Kante zusammengefasst, sodass das eine 
Flachenelement unmittelbar in das andere Flachenelement 
ubergeht. Da samtliche Flachenelemente eben sind, gehen 
sie nicht stetig ineinander uber und es entsteht ein 
facettierter Flachenverbund. 

Dieser facettierte Flachenverbund stellt ein Flachen- 
modell des digitalisierten Objekts dar. Wenn der Flachen- 
verbund der Oberflache ein endliches Volumen quasi 
wasserdicht umschliefit, so entsteht ein Volumenkorper- 
modell des digitalisierten Korpers. 

Das Flachen- oder Volumenk6rpermodell lasst sich 
problemlos durch CAD/FEM-Kopplung in ein FE-Programm 
importieren und zu einem FE -Model 1 vernetzen, mit dessen 
Hilfe physikalische Berechnungen durchgeftihrt werden 
kSnnen . 




Die Erfindung hat den Vorteil, dass sie fttr den Obergang 
zwischen der numerischen und analytischen Beschreibung 
der Objektdaten Kurvengleichungen ersten Grades heran- 
5 zieht. Im Vergleich zum Reverse Engineering, bei dem 
Kurvengleichungen dritten oder hoheren Grades verwendet 
werden, reduziert sich dadurch die Komplexit&t des zu 
losenden Gleichungssystems erheblich. Auch dann, wenn von 
einem sehr feinen Netz ausgegangen wird, urn ein besonders 

10 genaues dreidimensionales Modell des interessierenden 

Objekts zu erzeugen, ist die durch das weniger komplexe 
Gleichungssystem erzielte Zeitersparnis so grofi, dass 
sich der Rechenaufwand trotz der hohen Genauigkeit des 
Modells insgesamt verringert. Mit Hilfe der Erfindung 

15 lasst sich also mit relativ geringeia Rechenaufwand ein 
genaues dreidimensionales Flachen-, Volumenkorper- oder 
FE -Modell eines real bestehenden Objekts erzeugen. 

Die Erfindung kann sowohl in Form eines Verfahrens als 
20 auch in Form einer Vorrichtung oder eines Software- 
programmprodukts umgesetzt werden, 

Bei der Vorrichtung ist zu beachten, dass neben einer 
Datenverarbeitungseinrichtung, die die Daten- 

25 verarbeitungsschritte Einlesen des Netzmodells, Zerlegen 
des Netzmodells in bilineare Fiachenelemente, Vereinigen 
der bilinearen Flachenel entente zu einem Flachen- oder 
Volumenkorpermodell und gegebenenf alls Erstellen eines 
FE-Modells aus dem Flachen- oder Volumenkorpermodell 

30 ausfuhrt, auch eine Digitalisierungseinrichtung 

vorgesehen ist, mit der sich das Netzmodell des Objekts 
erzeugen lasst. Unter eine solche Digitalisierungs- 
einrichtung fallen sSmtliche bildgebenden Gerate wie 
Fotoapparate und Rontgengerate, die zweidimensionale 

35 analoge oder digitale Bilder erzeugen, aus denen in 
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Kombination mit einer Bildverarbeitung ein dreidimensio- 
nales Netz gewonnen werden kann. Unter eine solche 
Digitalisierungseinrichtung fallen aber auch optische und 
beruhrende Abtastgerate, die die Oberflache des Objekts 
abtasten, urn eine dreidimensionale Punktwolke zu 
erzeugen, aus der dann in Kombination mit einer Bild- 
verarbeitung ein Netz gewonnen wird. In diesem Zusammen- 
hang ist es ohne. Belang, ob die Digitalisierungs- 
einrichtung und die Datenverarbeitungseinrichtung 
raumlich getrennt sind oder ob die Digitalisierungs- 
einrichtung die Datenverarbeitungseinrichtung zur 
Ausfuhrung der Bildverarbeitungsschritte nutzt. 

Es liegt auf der Hand, dass die Erfindung anstatt durch 
da.e Vorrichtung auch durch ein Computerprogrammprodukt 
realisiert werden kann, das die angesprochenen Daten- 
verarbeitungsschritte anhand von Softwareroutinen 
abarbeitet, wenn es auf einem Computer lauft. Das 
Computerprogrammprodukt kann auf einem Datentrager 
gespeichert sein oder direkt in den Arbeitsspeicher des 
Computers geladen werden. 

Weitere Aufgaben, Vorteile und Merkmale der Erfindung 
ergeben sich aus der folgenden ausfuhrlichen Beschreibung 
von Ausfuhrungs- und Vergleichsbeispielen. Dabei wird auf 
die beigefugten Zeichnungen Bezug genommen, die Folgendes 
zeigen: 

Fig. 1 ein durch optische Abtastung digitalisiertes 
Abbild eines Probandengesichts in Form einer Punktwolke; 

Fig. 2 die Punktwolke von Fig. 1, nachdem sie ausgedtinnt 
wurde; 

Fig. 3 ein aus der Punktwolke von Fig. 2 erzeugtes Netz; 
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Fig. 4 ein Ausschnitt aus dem in Fig. 3 gezeigten Netz; 

Fig. 5a bis Fig. 5d eine erlauternde Darstellung, wie aus 
drei Polygonen des in Fig. 4 gezeigten Ausschnitts drei 
bilineare Flachenelemente gebildet land wieder vereinigt 
we r den; 

Fig. 6 ein aus dem in Fig. 3 gezeigten Netz erzeugtes 
Flachenmodell ; 

Fig. 7 ein aus dem in Fig. 6 gezeigten Flachenmodell 
erzeugtes FE -Model 1 der Gesichtsweichteile; 

Fig. 8 eine digitalisierte RSntgenschichtaufnahme eines 
menschlichen Schadels in H6he des Unterkiefers, in denen 
die Schadelgrenzen markiert und mit Punkten versehen 
sind; 

Fig. 9 ein aus mehreren Schichtaufnahmen erzeugtes 
digitales Abbild eines Schadels in Form einer Punktwolke, 
nachdem sie homogenisiert wurde; 

Fig. 10 ein aus der Punktwolke von Fig. 9 erzeugtes Netz; 

Fig. 11 ein aus dem in Fig. 10 gezeigten Netz erzeugtes 
Volumenkorpermodell des Schadels; 

Fig. 12a ein Netzmodell des menschlichen Ohres und 
Fig. 12b ein aus diesem Netzmodell durch Reverse 
Engineering erzeugtes Modell aus NURBS-Flachen; 

Fig. 13 ein durch Nachkonstruktion erzeugtes FE -Modell 
von Gesichtsweichteilen gemafi Stand der Technik; 
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Fig. 14 ein durch punktbasierte, schichtweise Vernetzung 
erzeugtes FE-Modeil eines Zahns gemafi Stand der Technik; 

Fig. 15 ein durch voxelbasierte, schichtweise Vernetzung 
erzeugtes FE -Model 1 eines menschlichen Schadels gemafi 
Stand der Technik; und 
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Fig. 16 ein durch Vernetzung einer Punktwolke erzeugtes 
FE -Model 1 eines Unterkiefers gemafi Stand der Technik. 

Anhand der Figuren 1-7 wird nun ein erstes Ausfiihrungs- 
bexspiel beschrieben, wie sich aus einem digitalisierten 
Abbild eines Probandengesichts ein dreidimensionales 
Flachenmodell und ein FE-Modell der menschlichen 
15 Gesichtsweichteile erzeugen lassen. 

Das Gesicht eines Probanden wurde mit Hilfe eines licht- 
codierten Triangulationsverfahrens (TRICOLITE™ der Firma 
Steinbichler) digitalisiert . Dabei wurde durch einen LCD- 
20 Pro D ektor eine Serie von Streifenmustern auf das Gesicht 
geworfen, die durch zwei CCD-Kameras aus unter- 
schiedlichen Blickwinkeln erfasst wurde. Der komplette 
Messvorgang dauerte etwa zwei Sekunden. Durch 
geometrische Auswertung (Triangulationsprinzip) wurde 
daraus ein dreidimensionales Abbild der Gesichts- 
oberflache in Form einer Punktwolke erzielt. Nahere 
Einzelheiten zu diesem lichtcodierten Triangulations- 
verfahren werden der Dissertation (citation in progress) 
von c. Holberg "Erfassung von Gesichtsoberf lachen durch 
ein lichtcodiertes Triangulationsverfahren", Ludwig- 
Maximilians-Universitat Munchen (2002) zu entnehmen .sein. 

Die erhaltene Punktwolke wurde anschliefiend gefiltert, urn 
eine bestimmte AuflSsung zu erzielen und redundante Daten 
emzusparen. Und zwar wurden die Bildpunkte, deren Lage 
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gegentiber den benachbarten Bildpunkten nur wenig abwich, 
geloscht, wahrend die Bildpunkte, deren Lage gegentiber 
den benachbarten Bildpunkten starker abwich, beibehalten 
wurden. Dadurch ergab sich die in Fig. 2 gezeigte 
ausgedtinnte Punktwolke, in der die Bildpunkte umso 
dichter liegen, je starker sich die Topologie der 
Gesichtsoberflache andert. 

Die ausgedtinnte Punktwolke wurde anschliefiend in das 
Bildbearbeitungsprogramm Rapid Form™ (INUS Technology, 
Inc.) importiert und zu einem Polygonnetz aus dreieckigen 
Polygonen vernetzt. Urn ein durchgehendes, sauberes 
Polygonnetz zu erhalten, wurde das Polygonnetz von 
LQchern und von kreuzenden, redundanten und nichtmannig- 
faltigen Flachen befreit. Im vorliegenden Fall wurde die 
Moglichkeit, das Polygonnetz zu homogenisieren oder die 
Anzahl der Polygonflachen zu erhohen oder zu reduzieren, 
nicht genutzt. Es entstand das in Fig. 3 gezeigte 
lochfreie und gereinigte Polygonnetz, das im DXF-Format 
zwischengespeichert wurde. 

Das im DXF-Format abgespeicherte Polygonnetz wurde 
anschliefiend in das Programm PolyTrans™ (Okino Computer 
Graphics) importiert, urn das Polygonnetz, ohne weiter 
verandert zu werden, im neutralen IGES-Format 
abzuspeichern, das den Datenaustausch zwischen unter- 
schiedlichen CAD-, CAM- und Computervisualisierungs- 
systemen erlaubt. Durch den Export als IGES-Datei wurde 
das Polygonnetz in bilineare NURBS-Patches mit der 
Elementnummer 128 zerlegt. Die entstandene IGES-Datei 
wurde anschliefiend in das CAD-Programm Mechanical 
Desktop- (Autodesk Inc.) eingelesen, wobei das Abbild des 
Probandengesichts nunmehr in Form von einzelnen Flachen- 
elementen (bilinearen NURBS-Patches) vorlag, die jeweils 
einem Polygon des ursprtinglichen Polygonnetzes 
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entsprachen. Mat der Funktion "Flachen zusammenhef ten" 
(in anderen Programmen auch als "Zusammenfiigen" oder 
"Stitchen" bezeichnet) wurden die einzelnen Flachen- 
elemente dann wieder zu einem Flachenverbund vereinigt, 
sodass das in Fig. 6 gezeigte Flachenmodell des 
Probandengesichts entstand. Dieses Flachenmodell wurde 
dann uber die CAD/FEM-Schnittstelle AMACISOXJT des 
Programms iia SAT-Format exportiert, urn die Geometrie des 
Flachenmodells der CAD/FEM-Kopplung zuganglich zu machen. 

Die einzelnen Verarbeitungsschritte, denen die Bilddaten 
beim Zerlegen des Polygonnetzes in die einzelnen Flachen- 
elemente und beim Vereinigen der Flachenelemente 
unterlagen, werden nun ausfuhrlicher anhand von Fig. 4 
15 und den Figuren 5a bis 5d erlautert. 

Fig. 4 zeigt einen Ausschnitt des in Fig. 3 gezeigten 
Polygonnetzes, der dem rechten Wangenbereich des 
Patienten entnommen wurde. In diesem Bereich werden nur 
die drei fett markierten Polygone betrachtet, deren Eck- 
bzw. Knotenpunkte hervorgehoben sind. Diese drei Polygone 
sind in Fig. 5a ohne ihre Umgebung dargestellt. 

Die in Fig. 5a dargestellten Polygone wurden durch die 
Konvertierung ins IGES-Format in die drei in Fig. 5b 
gezeigten bilinearen NURBS-Patches umgewandelt. Unter 
NURBS-Patches sind wie gesagt Flachenelemente zu 
verstehen, die jeweils durch zwei nicht-uni forme 
rationale B-Splines definiert sind, d.h. durch zwei frei 
formbare polynomiale Kurven. Da die Polygone, von denen 
ausgegangen wurde, eben waren, sind die beiden B-Splines, 
die jedes Flachenelement definieren, ebenfalls nicht 
gekrummt und entsprechen somit Kurven erstes Grades. Bei 
der Erfindung wird also jedes NURBS-Patch nicht mehr 
numerisch durch die Eckpunkte des jeweiligen Flachen- 
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elements, sondern analytisch durch zwei Kurven ersten 
Grades bzw. durch zwei Strecken beschrieben. Da die 
ursprtinglichen Polygone dreieckig waren, sind auch die 
entstandenen NURBS-Flachen dreieckig. Das bedeutet, dass 
5 jeweils ein Endpunkt der beiden Strecken des NURBS- 
Patches in einem Punkt zusammenf ailt und die beiden 
Strecken zwei Kanten des dreieckigen NURBS-Patches 
bilden, wahrend "sich die dritte Kante durch Verbinden der 
beiden verbliebenen Streckenendpunkte ergibt. Da jeder 
10 NURBS-Patch durch ein eignes Paar linearer B-Splines 
definiert ist, hat jedes NURBS-Patch seine eigenen 
Kanten, die es nicht mit den angrenzenden Fiachen- 
elementen teilt. 

15 Die einzelnen NURBS-Patches wurden dann durch die 

Zusammenhef tfunktion wie in Fig. 5c und Fig. 5d gezeigt 
jeweils mit den benachbarten NURBS-Patches vereinigt. Das 
Zusammenhef ten erfUllt dabei zwei Aufgaben: Zum einen 
werden mit dieser Funktion zwei oder mehr zusammen- 

20 hangende Flachen zusammengehef tet, urn einen Flachen- 

verbund zu erstellen, und zum anderen werden Fehler in 
der Geometrie oder Topologie, die wahrend der 
Konvertierung aufgrund von unterschiedlichen internen 
Toleranzen und Berechnungsf ehlern auftreten k6nnen, 

25 korrigiert. Durch das Zusammenhef ten entsteht also ein 
durchgangiger Flachenverbund, der als FlSchenmodell 
herangezogen werden kann. 

Das in Fig. 6 gezeigte, liber die CAD/FEM-Schnittstelle 
30 exportierte Fiachenmodell des Probandengesichts wurde in 
das FE-Programm Design Space™ von Ansys, Inc. importiert. 
Das Flachenmodell wurde dabei wie ein gekrtimmtes Fiachen- 
tragwerk behandelt, das sich nach der Schalentheorie 
verhait. Nach Zuweisung einer einheitlichen Dicke wurde 
35 das Konstrukt zu dem in Fig. 7 gezeigten dreidimensiona- 




len FE-Modell der Gesichtsweichteile vernetzt. Die 
Vernetzung erfolgte ohne Schwierigkeiten, da fehlerhafte 
Fiachen bereits wahrend der Reinigung des Polygonnetzes 
beseitigt wurden, sodass es zu keinen ttberschneidungen 
5 kam. Das FE-Modell wies im Grofien und Ganzen die gleiche 
hohe Auflosung wie das Polygonnetz und das im CAD- 
Programm erstellte Flachenmodell auf . Durch Festlegung 
entsprechender Lagerungen und Lasten liefien sich mit 
diesem FE-Modell hochauf losend Verf ormungen, Spannungen 
10 und Dehnungen in den Gesichtsweichteilen berechnen. Diese 
Berechnungen konnen zum Beispiel ftlr die Planung 
kosmetischer Operationen genutzt werden. 

Anhand der Figuren 8-11 wird nun ein zweites Ausfuhrungs- 
15 beispiel beschrieben, wie sich aus digitalisierten 

Schichtbildern eines menschlichen Schadels ein drei- 
dimensionales Voluiaenkorpermodell erzeugen lasst. 

Von dem Schadel eines Probanden wurde mit Hilfe eines 
20 Computertomograf ieverfahrens ein Satz von zweiundvierzig 
digitalen Ront gens chichtbil der n gewonnen. Mit Hilfe des 
Bildbearbeitungsprogramms 3D-Doctor' IM (Able Software) 
wurden in den Schichtaufnahmen die Grenzen des Schadel- 
knochens identif iziert und die markierten Grenzlinien mit 
25 mehreren Punkten versehen. Fig. 8 zeigt ein Beispiel 
eines solchen mit Grenzlinien und Punkten versehenen 
Schichtbilds im Bereich des Unterkief erknochens . 

Die in den einzelnen Schichtbildern markierten Punkte 
30 wurden dann zu einer dreidimensionalen Punktwolke 

zusammengef asst, die redundanten Daten durch Filtern mit 
einer bestimmten Aufl6sung entfernt und die verbleibende 
Punktwolke homogenisiert . Daraus ergab sich die in Fig* 9 
gezeigte ausgedtlnnte Punktwolke, in der die Bildpunkte 
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umso dichter liegen, je starker sich die Topologie der 
Schadeloberflache andert. 

Die ausgedilnnte Punktwolke wurde anschlieJJend zu einem 
Polygonnetz aus dreieckigen Polygonen vernetzt und das 
Polygonnetz von kreuzenden, redundanten und nichtmannig- 
faltigen Flachen gereinigt. Eventuell entstandene Locher 
wurden wieder geschlossen. Das sich ergebende Polygonnetz 
ist in Fig. 10 gezeigt. 

Das Polygonnetz wurde dann ahnlich wie im ersten 
Ausfuhrungsbeispiel in das Programm PolyTrans™ (Okino 
Computer Graphics) importiert und in bilineare NURBS- 
Patches mit der Elementnummer 128 zerlegt, die 
anschliefiend wieder unter Korrektur der aufgetretenen 
Geometrie- oder Topologiefehler zu einem Flachenverbund 
zusammengeheftet wurden. Da bereits das Polygonnetz frei 
von Lochern war und der Flachenverbund daher der 
durchgehenden Oberflache eines in sich geschlossenen 
festen Korpers (Solid) entsprach, wurde beim Zusammen- 
heften automatisch ein Volumenkerpermodell erzeugt. Das 
fertige Volumenkorpermodell ist in Fig. 11 gezeigt und 
hat im Grofien und Ganzen die gleiche hohe Aufldsung wie 
das ursprtingliche Polygonnetz, sodass es far ein hoch- 
aufiesendes FE -Mode 11 geeignet ist. Das entstandene FE- 
Modell konnte zur Simulation der Auswirkungen von 
Gewalteinwirkungen auf den Schadel herangezogen werden. 

Die obigen Ausfuhrungsbeispiele zeigen, dass sich mit der 
Erfindung hochgenaue dreidimensionale Modelle real 
bestehender Objekte erzeugen lassen. 

Ein Vergleich der Erfindung mit dem herkdmmlichen 
Verfahren des Reverse Engineering soli aufierdem den 



Nachweis bringen, dass sich mit der Erfindung auch der 
Rechenaufwand senken lasst. 

Hierzu wurde von einem aus 338 6 dreieckigen Polygonen 
bestehenden Netzmodell des menschlichen Ohres ausgegan- 
gen. Das in Fig. 12a gezeigte Netz wurde zum einen nach 
den Maflgaben der Erfindung und zum anderen mit Hilfe der 
Reverse-Engineering-Funktion des Programms Rapid Form™ 
(INUS Technology, Inc.) in NURBS-Patches der IGES- 
Elementenummer 128 umgewandelt. In beiden Fallen wurden 
die 3386 Polygone in ebenso viele NURBS-Patches 
umgewandelt. Im Fall des Reverse Engineering wurde 
aufierdem ein Modell mit einer geringeren Aufl6sung von 
232 NURBS-Patches angefertigt, das in Fig. 12b gezeigt 
ist. Die Berechnungszeit und der Speicherbedarf fur die 
einzelnen Modelle lassen sich Tabelle 1 entnehmen. 



Tabelle 1 





NURBS 
Anzahl 


Rechen- 
zeit 


Speicherbedarf 
(imkorcpriiriiert) 


Speicherbedarf 
(kartprimiert) 


Erfindung 


3386 


2,9 s 


1,2 MB 


0,4 MB 


Reverse 
Engineering 


3386 


1983 s 


29,6 MB 


16,1 MB 


Reverse 
Engineering 


232 


98 s 


2,0 MB 


1,1 MB 



Tabelle 1 ist zu entnehmen, dass das erf indungsgemafie 
Verfahren verglichen mit dem Reverse Engineering bei 
gleicher Auflosung (d.h. bei gleicher Anzahl an NURBS- 
Patches) deutlich weniger Rechenzeit und Speicherbedarf 
erfordert. Wenn beim Reverse Engineering die Anzahl an 
NURBS-Patches gesenkt wird, lassen sich zwar die Rechen- 
zeit und der Speicherbedarf reduzieren, doch werden 
dadurch geometrische Ungenauigkeiten in Kauf genpmmen. 
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Wie der Vergleich der Figuren 12a und 12b ergibt, erzeugt 
das Reverse Engineering zwar sehr glatte Oberflachen- 
ubergange, doch kann es insbesondere im Randbereich des 
Modells zu fehlerhaften Darstellungen kommen (siehe Pfeil 
in Fig. 12b) . Die Ungenauigkeiten im Randbereich sind 
umso auffalliger, je komplexer die nachzubildende 
Struktur ist und je geringer beim Reverse Engineering die 
AuflSsung (Anzahl der erzeugten NURBS-Patches) ist. 
Vergleichbare Darstellungsfehler treten bei der Erfindung 
nicht auf. 



Daruber hinaus stellte sich heraus, dass die durch das 
Reverse Engineering erzeugten NURBS-Patches insbesondere 
dann, wenn das Patch-Layout autoiaatisch erstellt wurde, 
zu irregularen, fehlerhaften und vertwisteten CAD-Flachen 
ftihrten, die nicht (z.B. uber CAD/FEM-Kopplung) weiter- 
verarbeitet werden konnten. Da die manuelle Erstellung 
des Patch-Layouts zeitaufwendig, fehlertrachtig und in 
der Regel nur vom Fachmann zu bewerkstelligen ist, bietet 
das Reverse Engineering keine praxis gerechte Alternative, 
urn aus digitalisierten anatomischen Strukturen wie dem in 
Fig. 12a gezeigten Ohr FE-Modelle zu erzeugen. 

Wie sich aus der obigen Beschreibung ergibt, eignet sich 
die Erfindung nicht nur als hochgenaues, Zeit und 
Speicherplatz einsparendes Alternatiwerfahren ftir das 
Reverse Engineering, sondern erschliefit die Erfindung 
auch neue Anwendungsbereiche bei der Modellierung 
hochkomplexer Objekte, etwa in den Biowisserischaf ten, die 
dem Reverse Engineering und anderen herkommlichen 
Verfahren bislang nicht zuganglich waren. 
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PATENTANSPRUCHE 

1. Verfahren zur Erzeugung eines dreidimensionalen 
Modells ftir ein real bestehendes Objekt, mit den 
Schritten: 

Digitalisieren des Objekts, urn ein Netzmodell des 
Objekts zu erzeugen; 

Zerlegen des Netzmodells in bilineare Fiachen- 
elemente; und 

Vereinigen der bilinearen Flachenelemente zu einem 
Flachen- oder Volumenkorpermodell des Objekts. 

2. Verfahren nach Anspruch 1, bei dem die bilinearen 
Flachenelemente dreieckig sind. 
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3. Verfahren nach Anspruch 1 oder 2, bei dem die 
bilinearen Flachenelemente NURBS-Patches sind. 



4. Verfahren nach Anspruch 1, mit dem Schritt Erstellen 
eines Finite-Elemente-Modells aus dem Flachen- oder 
Volumenkorpermodell . 
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5. Verfahren nach Anspruch 1, bei dem das Netzmodell 
aus Punktwolkendaten des Objekts gewonnen wird. 

5 6. Vorrichtung zur Erzeugung eines dreidimensionalen 
Modells fur ein real bestehendes Objekt, mit: 

einer Digitalisierungseinrichtung, mit der sich ein 
Netzmodell des Objekts erzeugen lasst; und 

einer Datenverarbeitungseinrichtung zur DurchfUhrung 
10 folgender Datenverarbeitungsschritte : Einlesen des 
Netzmodells; Zerlegen des Netzmodells in bilineare 
Flachenelemente; und Vereinigen der bilinearen Fiachen- 
elemente zu einem Fiachen- oder Volumenkdrpermodell des 
Objekts. 

15 

7 • Vorrichtung nach Anspruch 6, bei der die Daten- 
verarbeitungseinrichtung aus dem Fiachen- oder Volumen- 
kSrpermodell ein Finite-Elemente-Modell des Objekts 
erstellt . 

20 

8. Vorrichtung nach Anspruch 6 oder 7, bei der die 
Datenverarbeitungsschritte in der Datenverarbeitungs- 
einrichtung anhand von Sof twareroutinen abgearbeitet 
werden . 

25 

9. Computerprogrammprodukt zur Erzeugung eines 
dreidimensionalen Modells ftir ein real bestehendes 
Objekt, das anhand von Sof twareroutinen folgende Daten- 
verarbeitungsschritte abarbeitet, wenn es auf einem 

30 Computer lauft: Einlesen eines Netzmodells des Objekts; 
Zerlegen des Netzmodells in bilineare Flachenelemente; 
und Vereinigen der bilinearen Flachenelemente zu einem 
Flachen- oder Volumenkorpermodell des Objekts. 
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10. Computerprogrammprodukt nach Anspruch 9, das die 
bilinear en Flachenelemente erzeugt, indem es das 
Netzmodell des Objekts durch Konvertierung ins IGES- 
Format in bilineare NURBS-Patches zerlegt. 

11. Computerprogrammprodukt nach Anspruch 10 , bei dem 
das ins IGES-Format konvertierte Netzmodell 
ausschlieiilich IGES- Flachenelemente der Nummer 128 
enthalt . 

12. Computerprogrammprodukt nach Anspruch 9, das aus dem 
Flachen- oder Volumenkorpermodell durch CAD/ FEM-Kopp lung 
ein Finite-Elemente-Modell des Objekts erstellt. 
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ZUSAMMENFASSUNG 

(Fig. 7) 

Zur Erstellung eines dreidimensionalen Modells ftir ein 
real bestehendes Objekt wird das interessierende Objekt 
zunachst digitalisiert, um ein Netzmodell des Objekts zu 
erzeugen. Anschliefcend wird das Netzmodell in bilineare 
5 Flachenelemente zerlegt und diese schliefilich wieder zu 
einem Fiachen- Oder Volumenkorpermodell vereinigt. Aus 
dem Fiachen- oder Volumenkorpermodell lasst sich durch 
CAD/FEM-Kopplung ein Finite-Elemente-Modell erstellen. 
Die erzeugten Modelle sind nicht nur genau, sondern ihre 
10 Erstellung erfordert auch relativ wenig Rechenaufwand. 
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